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Abstract: One proposal for deriving effective cosmological models from theories of quan¬ 
tum gravity is to view the former as a mean-field (hydrodynamic) description of the latter, 
which describes a universe formed by a ‘condensate’ of quanta of geometry. This idea 
has been successfully applied within the setting of group field theory (GFT), a quantum 
field theory of ‘atoms of space’ which can form such a condensate. We further clarify the 
interpretation of this mean-field approximation, and show how it can be used to obtain a 
semiclassical description of the GFT, in which the mean field encodes a classical statistical 
distribution of geometric data. In this sense, GFT condensates are quantum homogeneous 
geometries that also contain statistical information about cosmological inhomogeneities. 
We show in the isotropic case how this information can be extracted from geometric GFT 
observables and mapped to quantities of observational interest. Basic uncertainty relations 
of (non-commutative) Fourier transforms imply that this statistical description can only 
be compatible with the observed near-homogeneity of the Universe if the typical length 
scale associated to the distribution is much larger than the fundamental ‘Planck’ scale. As 
an example of effective cosmological equations derived from the GFT dynamics, we then 
use a simple approximation in one class of GFT models to derive the ‘improved dynamics’ 
prescription of holonomy corrections in loop quantum cosmology. 

Keywords: Models of Quantum Gravity, Cosmology of Theories beyond the SM, Lattice 
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1 Introduction 

The most promising hope for connecting theories of quantum gravity with observation is 
to understand their consequences for early universe cosmology. Conversely, input from 
quantum gravity may be needed in order to complete or replace the standard paradigm of 
inflation. Open theoretical issues of inflation, such as the choice of initial conditions or the 
origin of the inflationary potential [1] or the need for a resolution of the past singularity of 
inflationary universes [2], should ultimately be addressed in a more complete framework. 
More generally, with a deeper understanding of quantum gravity, one might obtain a dif¬ 
ferent physical picture for the very early universe that does not require introducing scalar 
fields, but is at the same time compatible with observations such as made by Planck [3]. 

In order to obtain clear predictions from quantum gravity-inspired approaches to cos¬ 
mology, it is essential to understand their precise relation to an underlying quantum gravity 
theory. One such approach is loop quantum cosmology (LQC) [4-6], in which the classical 
Big Bang singularity is resolved [7] and the issue of initial conditions for inflation can be 
addressed [8], but whose relation to the full theory of loop quantum gravity (LQC) [9-11], 
or the closely related group field theory (GFT) [12, 13], a second quantised version of LQC 
[14], is not fully clear. In LQC, a symmetry reduction to a minisuperspace model is per¬ 
formed before quantising the remaining degrees of freedom with LQC techniques, leading 
to quantisation ambiguities and obscuring the potential embedding into the full theory. 

One perspective on addressing this issue, advocated in [15], is to view LQC and minisu¬ 
perspace quantum cosmology as ‘single-patch theories’ in which an elementary small chunk 
of space is quantised in order to capture the dynamics of an exactly homogeneous universe. 
One expects a more complete picture, rich enough to also capture inhomogeneities, to arise 
from a ‘many-patch theory’ in which many such chunks can interact. Depending on the 
physical length scales one associates with these chunks, one might think of this as related 
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to the separate-universe approach in cosmology [16]. This perspective then calls for an ap¬ 
plication of concepts and methods from many-body quantum systems in condensed matter 
physics to quantum cosmology, as it suggests thinking of a macroscopic universe as a con¬ 
densate of many such ‘atoms’ of space. This idea had been discussed in various contexts, 
from the perspective of analogue gravity [17] as well as in quantum gravity [18, 19]. It was 
then explored in [20] in a lattice model starting from classical general relativity. 

The viewpoint that a cosmological universe arises from the condensation of many 
‘atoms of space’ is most naturally investigated in the setting of group field theory (GFT), 
which provides a quantum held theory language for discrete (simplicial) geometry in which 
the concept of a condensate can be made sense of: a GFT condensate dehnes a non- 
perturbative ground state, describing a phase away from the Fock vacuum around which 
perturbative physics is dehned in LQG. This Fock vacuum is analogous to the Ashtekar- 
Lewandowski vacuum [21] of LQG, and the excitations around it form the vertices of LQG 
spin networks. This vacuum corresponds to a degenerate geometry (zero expectation value 
for areas, volumes etc), and a non-degenerate continuum must be sought away from it. 

The idea of deriving quantum cosmology models from the dynamics of a GFT conden¬ 
sate was explored in a series of papers [22-25] (see also the related work [26-29]). In [22, 23], 
approximating a physical state with a (generalised) coherent state of the GFT held opera¬ 
tor, and hence working, in the simplest case, in a mean-held approximation, it was shown 
that Schwinger-Dyson equations encoding the GFT dynamics reduce to nonlinear, nonlocal 
differential equations for the mean held, or ‘condensate wavefunction’. Furthermore, in the 
simplest approximation, one obtains a linear equation that resembles a Wheeler-DeWitt 
equation in standard quantum cosmology. In a WKB limit for the isotropic case, one can 
reproduce the Friedmann equation of classical vacuum GR. One might view these results as 
suggesting a direct derivation of quantum cosmology models from the dynamics of suitable 
GFT Fock states, reproducing exactly the usual formalism in terms of a ‘wavefunction of 
the universe’ annihilated by a Hamiltonian constraint operator [30]. 

However, it was clear that this most direct equivalence of the effective condensate 
dynamics and a minisuperspace quantum theory cannot be expected. The ‘condensate 
wavefunction’ is not an actual wavefunction; one can think of it as a classical held, whose 
amplitude and phase correspond, in the case of a Bose-Einstein condensate, to classical 
properties such as the density and velocity of the condensate [31]. Its dynamics represent 
the collective, hydrodynamic description of the condensate. The equations dehning these 
dynamics are nonlinear, both for real condensates {e.g., the Gross-Pitaevskii equation) 
and for GFT condensates, and so one would have to make sense of nonlinear quantum 
mechanics. It was already understood in [24, 25] that the WKB approximation for the 
resulting dynamics does not capture a physically meaningful limit, even though it formally 
corresponds to li —^ 0: in a Bose-Einstein condensate, in this limit one would consider a 
huid with almost constant density but very high velocity; for GFT condensates, one would 
assume that the individual ‘atoms’ are very large and behave semiclassically individually. 

In this paper, we further explore and clarify the interpretation of GFT condensates 
in cosmological terms as outlined in [25], using not a standard (‘first-quantised’) Wheeler- 
DeWitt equation or a quantum cosmology wavefunction but focussing instead on expecta- 
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tion values of Fock space operators that are meaningfully defined on any GFT many-particle 
state. These expectation values are given a cosmological interpretation which is used to ob¬ 
tain cosmological dynamics from dynamical relations between these expectation values. We 
arrive at an interpretation in which, while GFT condensates are homogeneous as quantum 
states, their effective classical (statistical) description in the hydrodynamic approxima¬ 
tion generally includes inhomogeneities. In this hydrodynamic picture, the single-particle 
wavefunction of the ‘ground state’ that all quanta are condensed into is reinterpreted as a 
statistical probability distribution. For this distribution to define an exactly homogeneous 
classical geometry, one would have to demand that condensation occurs into an extremely 
peaked state, analogous to a delta distribution in position or momentum space for a Bose- 
Einstein condensate.^ This requirement is an additional condition on top of condensation, 
which cannot be satisfied in general, as it depends on the details of the dynamics. 

In general, the mean field has a finite spread on minisuperspace, and thus the hy¬ 
drodynamic approximation leads to a statistical distribution of classical ‘patches’ with 
different geometric data. Global information about the inhomogeneities, corresponding to 
moments f d^x in the approximate continuum description, can then be extracted 

from ‘global’ expectation values of GFT operators. These moments are independent for 
different n, and knowing them for many different n allows (in principle) the reconstruction 
of the inhomogeneities ^jJ{x) (in the simplest case of only isotropic scalar perturbations)^. 
We give an example in which, in the homogeneous case, the improved dynamics prescrip¬ 
tion for holonomy corrections in LQG [32] can be obtained from such expectation values. 
In general, the resulting equations depend on the inhomogeneities. We also show that, as 
there is a fundamental limit on the minimal spread of the mean field in metric variables, 
a statistically almost homogeneous classical metric (i.e. a very sharply peaked statistical 
distribution on minisuperspace) can only arise if the average areas in this distribution are 
large compared to the fundamental ‘Planck’ scale set by the parameter k of the GFT. Ex¬ 
plaining why this is the case is then the restatement in GFT condensate cosmology of the 
puzzle of how to explain the smallness of inhomogeneities in our own Universe. This inter¬ 
pretation also explains the shortcomings of the analysis of [28] in which it was suggested to 
add inhomogeneities by hand and interpret the resulting ‘condensate wavefunction’ directly 
as a wavefunction for background and inhomogeneities d la Halliwell-Hawking [33]. The 
hydrodynamic approximation already contains statistical inhomogeneities. 

2 Cosmology with group field theory (GFT) condensates 

This section provides a self-contained overview of GFT, condensate states in GFT and 
their relation to quantum cosmology. For more details on the general GFT formalism see 
e.g. [12, 13]; full details of the construction of GFT condensates are given in [23]. 

^This sense of exact homogeneity can at best be satisfied for either metric or connection variables, due 
to uncertainty relations for canonically conjugate variables. For the purposes of this paper, we will focus 
on the notion of homogeneity in the reconstructed classical metric. 

^Concretely, given n such integrals, one chooses a suitable n-dimensional function space with basis {/n} 
and reconstructs the Un in ip{x) = ’^nfnix)', a possible choice is, e.g., fn{x) = sin(nfco|i:|)/(nfco|a;|). 
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Group field theory (GFT) provides a quantum field theory language for simplicial 
geometry and for the kinematics and dynamics of LQG. It can also be viewed as a general¬ 
isation and extension of matrix models which achieve a definition of 2d quantum gravity in 
terms of random matrices [34]. In the formulation used so far in the construction of GFT 
condensates, one uses a complex scalar field 

ip-.G^/G^C ( 2 . 1 ) 

where G is a Lie group fixed from the outset, which will become the gauge group of gravity. 
For models of M quantum gravity, the conventional choices are G = SL(2, C), G = Spin(4) 
or G = SU(2), the latter being the gauge group of the Ashtekar-Barbero formulation of 
gravity underlying LQG. In this paper we will often assume G = SU(2). 

<y9 is a function on G^/G where G acts by diagonal right action, i.e. a function on four 
copies of G satisfying 

= 'ih^G. ( 2 . 2 ) 

In the geometric interpretation of an elementary GFT quantum as a tetrahedron, gj de¬ 
fine parallel transports along four links through the tetrahedron’s faces and (2.2) is the 
invariance under a gauge transformation acting on the vertex where the links meet. 

The dynamics is then defined by a choice of action S[ip,(p] and the functional integral 

/ n 

Vf, ^ Z[r] (2.3) 

where Aj are the coupling constants of the theory and the second equality is the perturbative 
expansion of Z (around (^ = 0) in Feynman graphs F, where each F is associated with a 
symmetry factor and a Feynman amplitude Z\r]. Each F forms a two-complex that can 
be interpreted as a discrete spacetime (with no boundary), or a spin foam. The amplitude 
Z[F] is then the spin foam amplitude defining the dynamics of LQG [11, 35]; there is a 
one-to-one correspondence (within a certain class of models) between spin foam models, 
defined by a choice of Z\r], and GFT models, defined by a choice of action S [36]. 

In addition to the covariant formalism for GFT in terms of the functional integral, 
there is also a canonical formalism that is closer to canonical LQG [14]. One defines a Fock 
space by starting with a Fock vacuum |0) annihilated by the field operator, (p{gi)ftf) = 0, 
and imposing canonical commutation relations 

W{gi),^{g'i)] = W\gi),^H9'i)] = O, [ip{gi),(p\gi)] = l{gi,gi) (2.4) 

as in usual non-relativistic bosonic field theory. \{gi,g'j) is an identity element on the space 
of fields compatible with (2.2); for compact G 

r ^ 

Hgi,gi) ■= / dL YiSig'jhgj^) (2.5) 

1=1 

where here and in the following dh is the normalised Haar measure on G. A one-particle 
state \gi) := ip^gi)\9) is then identified with a geometric tetrahedron, or an open 4- 
valent spin network vertex in LQG; more complicated simplicial geometries or LQG spin 
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networks correspond to many-particle Fock states. For an explicit construction of this 
Fock space starting from LQG spin-network wavefunctions, in complete analogy to the 
usual introduction of second quantisation starting from A^-particle wavefunctions, see [14]. 

One can use this Fock space construction to define condensates in analogy with con¬ 
densed matter physics. Such condensates can be characterised by a non-zero expectation 
value for the field operator, the GFT field and a large number of particles (potentially 

infinite) with respect to the Fock vacuum |0). They are usually defined as superpositions 
of states with different particle number. In the context of LQG, this means that one is not 
working on a fixed graph or discretisation of space, but taking a superposition of many 
discrete geometries into account. 

The simplest states used to describe such a condensate are of the form 

jcr) := AA(cr) exp (d^)'^ a{gi) (p\gi)^ |0), (2.6) 

which is a coherent state, an eigenstate of the field operator (p{gi) with eigenvalue (7{gi). 
As discussed in more detail below, the state (2.6) describes a mean-field approximation 
where the (non-zero) mean field is given by the ‘condensate wavefunction’ cr(gj). 

One now observes that the domain space G^/G is a finite-dimensional space of geomet¬ 
ric configurations, here the possible configurations of a tetrahedron, given in terms of the 
parallel transports of a (discrete) connection. This space still contains gauge-variant data, 
as one has to take into account gauge transformations acting from the left (on the open 
ends of the spin network links). In LQG, the gauge-invariant state is defined on a ‘dipole’, 
a graph given by four links and two vertices such that all links connect the two vertices, and 
the space of wavefunctions is L^(SU(2)\SU(2)^/SU(2)). Accordingly, in (2.6) we require 
left-invariance of the mean field cr(g/), cr{gi) = a{hgi) {a is already right-invariant because 
of (2.2)). For G = SU(2), the space G\G‘^/G is six-dimensional; a parametrisation of 
the resulting 12-dimensional phase space in terms of geometric variables corresponding to 
metric and connection is given in [24]. Viewing the tetrahedron as a locally homogeneous 
patch of the universe, one can view this phase space as the space of homogeneous (intrin¬ 
sic and extrinsic) geometric data, i.e. of homogeneous 3-metrics and Ashtekar-Barbero 
connections or second fundamental forms. 

For the GFT formalism [37], as for LQG [38], there exists a dual ‘flux representation’ 
which is obtained using a non-commutative Fourier transformation on G defined by [39] 

^{Bi) ■■= J{dg)^ ^neg,(B/)^ ip{gi) (2.7) 

where Bj are elements of the Lie algebra of G, and eg{B) is a choice of plane waves 
on G, satisfying standard properties such as compatibility with the adjoint action of G, 
eg{hBh~^) = e/j-ig/j(il), and completeness f dB eg(B) = 6{g). A standard choice for SU(2) 
is eg{B) = exp(-^tT(gB) / k) (where k and B have physical dimensions of area). In general, 
the choice of plane waves corresponds to a (necessary) choice of quantisation map for the 
fluxes [40], and the flux representation depends on this choice. Functions in the image of 
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the Fourier transform are then equipped with a non-commutative ★-product defined by 


eg[B) ★ eg'(B) := egg'(B) (2.8) 

on plane waves and extended by linearity. In particular, from (2.2) one hnds that 

‘f{Bi) = J{dg)^dh ^ eg^(.B/) ★ e/i(S/)^ ip{gj) = ip{Bj) Bj^ (2.9) 

with the non-commutative delta function (5*(i?) := f dg eg(B); the field (p(Bf) satisfies the 
(non-commutative) closure constraint Bj = 0. If the Bj are identified with bivectors 
representing integrals of an area 2-form over the faces, Bf^ Iaj e^ A e^, the closure 
constraint means that the faces close to form a tetrahedron. The dual Lie algebra vari¬ 
ables (corresponding to ‘fluxes’ of the triad in LQG) hence represent a (discrete) metric, 
canonically conjugate to the discrete connection given by g/. 

The connection to cosmology is made by rewriting the GFT quantum dynamics as 
dynamical equations for the mean held o'(g/) or its non-commutative Fourier transform 
d{Bi) and by constructing macroscopic geometric observables out of the elementary group 
and Lie algebra variables. Such observables are identihed with cosmological variables; for 
instance, the total volume in a GFT condensate dehnes the volume of a region of the 
universe, if a is the scale factor and Vq the coordinate volume. Such a region can be 
thought of as the ‘hducial cell’ in LQG [4-6] (which, as in LQG, can make up the whole 
universe for a compact universe), and itself consists of a large number of fundamental GFT 
quanta. Constraints between the cosmological variables dehne the effective cosmological 
description of the GFT dynamics (an effective ‘Friedmann equation’). 

In this paper we make this correspondence between GFT observables and cosmological 
variables, the crucial step in the physical interpretation of the effective GFT dynamics, more 
precise. We revisit the previous constructions, mainly given in [23, 25], and argue that, in 
the hydrodynamic approximation, the effective continuum classical geometry reconstructed 
from a GFT condensate is not necessarily homogeneous. We show how different GFT 
observables can be used to extract statistical information both about the homogeneous 
mode and the inhomogeneities. These points are general, and we illustrate them only in 
the simplest case of an isotropic universe where they are clearest. This simplest case will 
already be sufficient to make two points: the apparent near-homogeneity of the observed 
Universe requires a statistical distribution of ‘patches’ that are ‘large’ on average; and for a 
homogeneous, isotropic universe we give, under a few simplifying assumptions, a derivation 
of the LQG ‘improved dynamics’ prescription for holonomy corrections. We will see how 
the presence of inhomogeneities affects the effective cosmological equations. 


3 Effective classical geometry from global GFT observables 

Part of the geometric information contained in any GFT Fock state, e.g. a condensate 
state of the form (2.6), can be expressed in terms of expectation values of suitable second- 
quantised operators. We focus on ‘one-body operators’ of the general form 

O := y (dg)'^ (dc/')^ <pHgi) 0{gi, g'l) (p{g'i). (3.1) 
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Such operators are extensions of operators 6 of the first-quantised formulation in that the 
matrix elements of d are inserted into the integral (3.1), 0{gj,g'j) := {gj\d\g'j) for single¬ 
particle states \gi) and \g'j). This is the standard procedure in second quantisation; for 
further discussion of GFT as a second quantised formulation of LQG, see [14]. 

Equivalently, given a representation of the matrix elements of the Hrst-quantised op¬ 
erator d in terms of single-particle wavefunctions, 

(T|d|T') = I(dg)^ T(^(DT')(5/) (3.2) 

where D is in general a differential operator on G‘^, (3.1) can be written as 

O := y {dg)"^ (p^gi){D(p){gi). (3.3) 


The two expressions are equivalent as jT) = f (dg)'^ in terms of the basis 

For example, for G = SU(2) and interpreting GFT Fock states as LQG spin networks, 
the one-body operators corresponding to ‘total fluxes’ are defined as 


b] ■= iK j (dfi')'^ (exp (r)t) gj) 


(3.4) 


where k is an area, the ‘Planck’ area (in kinematical LQG, k = SvryhGN with the Barbero- 
Immirzi parameter 7 and Newton’s constant Gn), and r* is a basis of the Lie algebra of 
SU(2), usually taken as r* = in terms of the Pauli matrices a*. 

The total area associated to all I-th. faces in a given state is the expectation value of 


Ai := K j{dg)‘^ (p^{gj)y^-Ag^ ^ (gj) 


(3.5) 


where Agj denotes the Laplace-Beltrami operator on SU(2), acting on the I-th argument 
of If. The square root is to be defined in terms of eigenvalues of —Ag^, in the sense of 
Dirac [41]. These eigenvalues are of the form —j{j + 1) for integer or half-integer j, giving 
the celebrated discrete area spectrum of LQG. More concretely, introducing the coordinate 
system on (one half of) SU(2) 

g = y/l — TT 1 — i(T • TT , |vf| < 1 (3.6) 


(which associates to 5 G SU(2) a Lie algebra element -nig] := \a ■ TT{g\ G su(2)), we have 
^/ = y / (^--^1 -7f|af + ^p{ni) , 

i/ = ^ y - 7r}7r})a/aj + 37r}5/ (^(tt/) , (3.7) 

where there is summation over repeated indices i,j, k but not over I which is fixed. If we 
identify the differential operator appearing in (3.7) with a first-quantised flux operator b\ 
as in (3.3), i.e. 

b\ =■ j{dg)^ ^Hgi){b]p){gi), (3.8) 
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it is easy to check that indeed 



Operators of the form (3.1) or (3.3) are defined on the whole Fock space, so that quan¬ 
tities such as ‘the expectation value of the total area’ can be associated to any state. For 
states of cosmological interest, e.g. a condensate defined by (2.6), such expectation values 
are identified with cosmological observables such as the scale factor, and relations between 
expectation values, derived from the fundamental quantum dynamics, become relations 
between cosmological observables. This allows the derivation of generalised Friedmann 
equations from the fundamental GFT dynamics, for specific choices of dynamics [25]. 

The analysis of [25] focussed on the simplest one-body operators on the GFT Fock 
space, the total flux (3.4) and a ‘total group coordinate’ 

ft/ := J (dg)"^ ^gi] ^Hgj)ip{gj) (3.10) 

which are the analogue of total momentum and ‘total position’ in the GFT context (the 
latter is in contrast to the physically more meaningful centre-of-mass position [25]). One 
could, however, also consider other operators, e.g. the areas (3.5), or 

«/ := J (dfi-)^ ^Hsi) ^ V>) {gj) ■ (3.11) 

a/ corresponds to the sum of the squared areas associated to each I-th face in a general 
many-particle state. This is a very different quantity from the square of the area (3.5), 


ft/ = y(dfi-)^ {dg')"^ <fH9l) 



i9l)'^H9'i) 



(3.12) 


which is not a one-body operator of the form (3.1). In general, one cannot expect a simple 
relation between expectation values of dj and Aj, although in the case of a particularly 
simple state such as (2.6) there is a relation, given that all expectation values are expressible 
in terms of integrals involving the mean field a. 

For the state (2.6), the expectation values of Af and d/ are given by 



(d/) = / (d^)^ a{gi) ( b]b]a \ {gi) 


(3.13) 

If cr(gj) was a quantum-mechanical wavefunction, these expectation values would corre¬ 
spond to expectation values (^) and {A^) for an operator A, and could be used to compute 
the variance {A^) — {A)"^ which contains information about the statistics of the observable 
A. But this is not the right interpretation in second quantisation, d/ is not A?^ and dehnes 
an independent observable for general many-particle states. (y{gi) is a mean field on the 
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‘minisuperspace’ parametrised by gj, describing the collective properties of the coherent 
state (2.6). \a{gi)\'^ can be interpreted as a classical number density of classical ‘patches’; 

ixc) ■= j Xc{gi)^K9l)0{9l'^ = k(5/)P (3.14) 

gives the expectation value for the number of quanta for which {< 7 /} are in C C {xc 
is the characteristic function of C). Note again the difference between hrst and second 
quantisation; the analogue of xc in quantum mechanics would correspond to a projective 
measurement with eigenvalues 0 and 1, whereas the set of eigenvalues of xc is No- 

In this interpretation, the mean field (7{gi) is used to give a classical statistical de¬ 
scription of the hydrodynamic approximation, in which a generic condensate state does 
not describe a perfectly homogeneous universe, but rather a distribution of patches with 
different values for geometric quantities such as e.g. curvature invariants constructed from 
the gj. An approximately homogeneous universe arises from a sharply peaked mean held. 

Let us make this more precise. In [25], for the simplest case of an isotropic condensate, 
the identihcation of GFT expectation values with cosmological observables was 

(A) ^H'^[9j])-^^{'^[9j])^ =■ Ti o? , 

(ft/) = (dfi()^ T:[gi\ ‘fHgj)<f{gj)^ =: NVj siii{lo uj) , (3.15) 

where Tj and Vj are dimensionless su(2) ~ elements of order one, a is the cosmological 
scale factor, w the spin connection, N := (N) the expectation value of the number operator 

^ j0H9 i)0{9i) (3.16) 

and lo is a dimensionless number corresponding to a choice of coordinate units. The overall 
scaling of 11/ with N implements the observation [25] that holonomies should be given by 
intensive observables and the factor inside the argument corresponds to a choice of 

coordinate system in which the condensate as a whole is extended over a hxed coordinate 
volume, which essentially corresponds to a hxed hducial volume in LQC. We can introduce 
a hducial volume Vq as in LQC and in the following write Iq = Vq . 

One issue with using (3.15) as cosmological variables is that one might expect Tj and 
Vj to be zero by SU(2) symmetry^ and hence prefer SU(2) invariant quantities such as 
the area operator Aj. But there is another issue with the consistency of the low-curvature 
limit: as // and Tlj/N are canonically conjugate at least approximately as j/fj <C 1, they 
should correspond to classical observables with Poisson brackets 

|sin((V[)/iV)^/^a;),a^| oc -^ + 0{oj). (3-17) 

But in the limit of small ui the left-hand side is (Vo/N)^/^{a;, a^}, and hence {w,a^} must 
scale with the number of quanta as {N/Vq)^^^, so that the resulting classical Poisson brack¬ 
ets depend on N and hence indirectly on quantities like the scale factor'^. This seems to be in 
^Thanks to Lorenzo Sindoni for pointing this out. 

“^We are excluding the possibility that is a universal constant independent of cosmological observables, 
which would seem hard to justify physically. 
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conflict with the Poisson brackets of classical general relativity for which {w, a^} oc 

[4]. In order to obtain Poisson brackets consistent with the classical limit at low curvature, 

the total area must be defined with a nontrivial scaling and we dehne 

i^i) = ('^5')'^ -■ (3.18) 

with the same Vq as above (the total coordinate volume associated to the condensate), so 
that is a physical area. Here I is hxed so that we focus on only one of the four 

areas, as should be sufficient for isotropic universes. For a heuristic picture giving some 
additional justification to the factor in (3.18), one can think of a large ‘fiducial cell’ 
composed of many elementary cells; one would associate the cosmological area to 

the area of one of the sides of this cell. However, the total area {Aj) does not give just the 
area of one of these sides but, since one is summing over all tetrahedra, overcounts by a 
factor The scaling is also consistent with constructions such as [43] rooted in LQC. 

While (3.18) is simply a dehnition of a, different from the previous one in (3.15), also 
assuming that the expectation value {Aj) defines an extensive observable [i.e. (Aj) oc N) 
seems to imply that oc N. We will conhrm this argument explicitly in section 5, and 
derive this relation from an approximation of the GFT dynamics. 

We can now introduce additional geometric observables for GFT condensates and, 
in the hydrodynamic approximation, identify these with cosmological (metric) variables. 
This identification will make the effective classical statistical distribution of inhomogeneities 
explicit. For simplicity, we provide this identification for the case of an isotropic condensate 
(given e.g. by (2.6) with a choice of (T{gi) that incorporates isotropy). These observables 
we consider are again expectation values of one-body operators on the GFT Fock space. 

The hydrodynamic approximation to the quantum dynamics replaces the full quantum 
properties of the state by a finite number of expectation values; it can be viewed as a 
replacement of the quantum state by a classical statistical distribution of classical particles 
with properties such as momentum and position in condensed matter physics and areas, 
angles, etc. in the case of GFT. For instance, in this semiclassical description the condensate 
has a given number of quanta S> 1, given by the expectation value of the number operator 
N, whereas the quantum state is really a superposition of states with different numbers of 
particles. We will further clarify this interpretation in section 4; for a general discussion 
of this statistical nature of the hydrodynamic approximation, see also [42]. Here we use 
the statistical distribution on minisuperspace, given by the mean field, to define further 
observables that allow us to extract information about cosmological inhomogeneities. 

In this effective classical picture obtained from taking expectation values, the mean 
field defines a statistical distribution of microscopic geometries for N classical patches. 
For large N, we can represent this continuous distribution approximately by a collection 
of N patches, each labelled by an index i and associated with an area H)-, such that 
the statistical distribution of the different values H)- corresponds to the distribution given 
by the mean field. One can think of the as each representing a sufficiently small 
region in minisuperspace such that the average number of patches in this region is one. In 
this approximation, the statistical distribution for a classical random variable Aj is well 
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represented by a single ‘universe’ made up of N patches with definite Within this 
approximation, (3.18) becomes 

{Ai) = a2ivV3y^2/3 = ^ 4 =: ^ (l - 2V’‘) , (3.19) 


where in the last equality we have defined a quantity in terms of A\. To avoid over¬ 
counting, the ‘perturbations’ are constrained to satisfy 'ip^ = 0, so that the scale factor a 
is obtained from the average over the A^ The factor {Vq/N)'^/'^ in (3.19) is the coordinate 
area associated with each patch. 

(3.19) is simply a dehnition of a geometric quantity ip^. For small ^p^ and isotropy, so 
that metric fluctuations can be captured by a single (volume) variable, one would expect 
to recover the usual formalism of linear cosmological perturbations [45] from this setting, 
where one would really think oi ^p as a gauge-invariant (Bardeen) potential. At this kine- 
matical level, the ip^ are not necessarily small in any sense. The precise connection between 
this GFT formalism and the setting of linear perturbations in cosmology, and the physical 
interpretation of a scale factor a obtained from such an average, will presumably only be 
clear for condensates for which ■0 <C 1 and one can treat tp as a linear perturbation. 

Again, in this approximation in which the statistical distribution over geometric data 
is approximated by N classical quantities Aj, the expectation value of the operator d/ 
defined in (3.11) is 


(«/) = 







(3.20) 


as = 0. We can now extract information about the inhomogeneities from the expec¬ 

tation values (A/) and (d/): for a large number of quanta well approximating a continuum. 



N ) 4 (4)2 



(3.21) 


which is expressible only in terms of expectation values of the condensate and Vq which 
defines a choice of coordinate units. 

Similarly, defining an operator 4 (with n > 3) by 



we can identify its expectation value with 


(3.22) 


iM) = 


E 


n) 


2n/3 


^2n 


1-2V^M = 


{Mr 

Nn-l 


+ 


Nj 


2n/3 


n 

m=2 


n 


m 


E 


(3.23) 

and hence (in the continuum approximation) with a weighted sum of all moments of the 
perturbation ip, of the form f d^x ^p{x)”^, for 2 < m < n. These moments all give indepen¬ 
dent statistical information on the function ip. In cosmology, f d^x '0(x)^ corresponds to 
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the total power spectrum, whereas the higher moments correspond to the total bispectrum, 
trispectrum etc., of the scalar perturbations given by i/j. These quantities are directly re¬ 
lated to cosmological observations (again, in the regime where inhomogeneities are small) 
which are also statistical in nature. The hydrodynamic approximation in GFT provides us 
with a classical statistical description of what is really a coherent, homogeneous quantum 
state of geometry. This idea bears striking resemblance to the mechanism in inflation [44] 
of a transition from a quantum state of the inflaton to a classical statistical description 
for inhomogeneities, which is then observed in the cosmic microwave background. Viewing 
quantum cosmology as the hydrodynamics of quantum gravity suggests a similar mecha¬ 
nism for quantum geometry [42]. 

4 Homogeneity and classicality conditions for GFT condensates 

In section 3 we have argued that an effective classical picture constructed from the hy¬ 
drodynamic approximation to a generic GFT condensate is to be interpreted as a classical 
inhomogeneous universe: the mean field cr{gj) contains statistical information both about 
the homogeneous mode (corresponding, in the isotropic case, to the scale factor a) and 
about inhomogeneities. In this section, in order to connect with previous work, we recall 
the discussion of [23, sec. 3] which aimed at constructing states that are candidates for 
semiclassical, exactly homogeneous universes, and see where it needs to be extended. 

The classical phase space variables of a single GFT quantum, interpreted as a geometric 
tetrahedron, specify a discrete metric and a discrete connection. The momenta conjugate 
to the group elements gi are four Lie algebra elements (bivectors) Bj, subject to a closure 
constraint Bj = 0, so that only three are independent. Identifying 

Bf^ = e^>^efei (i = 1,2,3) (4.1) 

defines a discrete triad and a discrete 3-metric by gij = efejA (where indices are 
contracted with the appropriate G-invariant tensor). Equivalently, one can define 

9ij '■= ^iY{^BiB2B^) BkmBin , Bij := B^ Bjab ■ ( 4 - 2 ) 

The gij parametrise the space of gauge-invariant discrete metric data for one tetrahedron. 

In the construction of [23], one considers an embedding of V S> 1 tetrahedra into 
a manifold in which the coefficients gij {i and j label the three edges of a tetrahedron 
emanating from the same vertex) specify a continuum metric expressed in a fixed frame 
and evaluated at this vertex, i.e. 


gij =: g{ei,ej){x) (4.3) 

where {ej} is the given frame and x the position of the given vertex in the embedding. Put 
differently, in the embedding the edge i is aligned with the vector field e*. The construction 
assumes that the continuum metric to be reconstructed is almost constant on the scale of 
the tetrahedra, and requires a choice of vector fields e*. The latter are fixed by using an 
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embedding into a manifold with topology Ai — &/X and a transitive group action by the 
(arbitrary but fixed) Lie group ©. This group action provides a natural choice of {oj}: 
take {e,} to be a basis of left-invariant vector fields, which is unique up to a choice of scale 
in the Lie algebra of © and a global 0(3) rotation. 

It is then clear that a large number of tetrahedra all with the same gij approximate a 
spatially homogeneous metric, i.e. a metric compatible with the group action by ©. 

As a conclusion from this geometric interpretation at the classical level, the following 
criteria for GFT states to describe spatially homogeneous classical geometries, at the scale 
of the GFT quanta, are given in [23]: first, the quantum analogue of the classical property 
of homogeneity, i.e. having the same g^j for all tetrahedra, is taken to be ‘quantum 
homogeneity’, the condensation of many quanta into the same microscopic quantum state, 
specified by the mean field cr{gi) in a coherent state such as (2.6)^; second, a semiclassicality 
condition is required for the ‘wavefunction’ cr{gj), later taken to be the validity of the 
WKB approximation when applied to (j{gi). Two further conditions — near-flatness of 
the geometry on the scale of the tetrahedra, and a large number of quanta N ^ 1 — are 
required for consistency for the reconstruction procedure and geometric interpretation. 

In order to elaborate on the interpretation of these criteria, it may be helpful to recall 
the physical meaning of the mean-field approximation and the mean field (y{gi) in the 
context of the physics of Bose-Einstein condensates. This is standard textbook material 
(and we follow the discussion of [31, ch. 2] closely), which may however be less familiar to 
practitioners of quantum gravity. 

In second quantisation in non-relativistic quantum mechanics, a field operator 'L(r) 
is introduced as a superposition of annihilation operators associated to a complete set of 
first-quantised wavefunctions, 

^(0 = X] ho + ^ a„ , (4.4) 

i' 

where i/ is a set of labels characterising the states and “0” denotes the ground state. In 
the simplest case of a non-interacting Bose gas in a box, the wavefunction (/>o is simply a 
constant and the higher (pi, are plane waves. In the Bogoliubov approximation, one then 
replaces the operator oq by the c-number y/No, where Nq is the number of atoms condensed 
into the ground state, and treats the rest as a small fluctuation, 

<po{f^ ao ^Q{r)-.=- s/Wq (po{r), 6i!{r) ■.= '^(py{r) hy . (4.5) 

The mean field approximation is then the limit in which the fluctuations hT are ignored. 
In the words of [31], in this approximation “[...] the field operator coincides exactly with 
the classical field Tq and the system behaves like a classical object. This is the analogue 
of the classical limit of quantum electrodynamics where the classical electromagnetic field 
entirely replaces the microscopic description of photons.” 

®This quantum notion of homogeneity, which is crucial to the condensate assumption but must be dis¬ 
tinguished from (statistical) homogeneity of the reconstructed classical geometry, is generalised to ‘wave- 
function homogeneity’ in the construction of generalised condensate states in [29]. 
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The A^o atoms condensed into the ground state given by 0o (^) are not necessarily found 
at the same point tq, or even close to some tq. For the non-interacting gas, they are just 
evenly distributed over the box. In this case, in momentum space <^o(^ oc 5{p) so that 
there is condensation into a single value (zero) for the momentum. This is true for the free 
gas; for a general interacting system, the ground-state wavefunction will have some finite 
spread both in momentum and position space, subject to the uncertainty relation 

^xAp > I. (4.6) 

The mean field \ho is not semiclassical in the WKB sense, having constant phase in the 
simplest case and slowly varying phase more generally. 

Coming back to condensates in GFT, we have implemented the mean-field approx¬ 
imation by an appropriate choice of state in the GFT Fock space. The simplest choice 
( 2 . 6 ) is simply an eigenstate of the field operator (p{gi) with eigenvalue cr{gi) so that under 
expectation values (with normal ordering) (p{gi) (^{gj) as in (4.5). This is then already 
the semiclassical description (hydrodynamic approximation) of the full GFT dynamics; the 
operator (p{gi) is replaced by the classical field (j{gi). 

On the other hand, requiring that the classical metric geometry one can reconstruct 
from the hydrodynamic approximation is exactly homogeneous, at the scale of the GFT 
quanta, is analogous to requiring the density of the fluid describing a Bose-Einstein con¬ 
densate in the hydrodynamic approximation to be peaked around a single point (in position 
space or momentum space). Here we would require that the analogue of the ground-state 
wavefunction, the mean field a when the GFT dynamics is imposed, is sharply peaked 
around one particular metric geometry; more concretely, the non-commutative Fourier 
transform d{Bi) of (j{gi) has to be sharply peaked around values of Bj that correspond, 
by (4.2), to the same 3-metric gij. This is consistent with our analysis in section 3: the 
magnitude of inhomogeneities depends on the shape of the function a in minisuperspace. 

We can then say that the assumption of condensation is a quantum notion of homo¬ 
geneity, and the hydrodynamic description arising from the mean-held approximation is a 
semiclassical approximation®. No further assumption of semiclassicality should be imposed 
on the mean held, as was already discussed previously in [24, 25]. 

Whether there is an approximate classical, homogeneous description of the condensate, 
in terms of the statistical distribution of patches reconstructed from the hydrodynamic 
approximation, depends however on the properties of the mean held. 

In order to make more precise statements about the magnitude of inhomogeneities, it 
is necessary to solve the (approximate) GFT dynamics to get a physical cr{gj). There is 
a fundamental limit on the minimum spread of a in the group and Lie algebra variables. 
Using the coordinate system on SU(2) as in (3.6), 

g = Vl — if 1 — ia ■ TT, jifl < 1, 

^Generally, outside of the mean-field approximation for coherent states, one can think of other ways of 
implementing different notions of either homogeneity or semiclassicality, or both. One example is given by 
the generalised condensates of [29], which are ‘quantum homogeneous’ in the sense of being determined by 
a single wavefunction for all quanta, but also involve coarse graining of microscopic degrees of freedom. 
Their hydrodynamic description is not simply a classical limit for the GFT quantum field. 
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and defining the non-commutative Fourier transform by 

/» := I d 5 = J f{n) (4.7) 

|7f|<l 

where B =: \a ■ B for B G su(2), we have the relation 

AvrAi? > I (4.8) 

where we use the requirement on GFT condensates to describe geometries that are nearly 
flat on the scale of the tetrahedra, i.e. to be peaked on values |7f| <C 1 (where vf are 
coordinates on a suitable gauge-invariant combination of group elements such as gig^^, 
see [24]), where (4.7) becomes the standard Fourier transformation. In general, for any 
/, Avr < 1 which gives a lower bound on AB] the image of the non-commutative Fourier 
transform consists of functions on su( 2 ) ~ with hnite resolution of order k. 

The uncertainty relation (4.8) has an interesting consequence. The magnitude of the 
inhomogeneities is controlled by the relative spread {AB)/{B). In order to achieve 
if} ~ I0“'^'® as required observationally, together with near-flatness Avr <C 1, we then need 
to have (B) ^ 10^'^k: the average (physical) length scale associated to the statistical distri¬ 
bution over geometries needs to be at least a few orders of magnitude larger than the scale 
set by the ‘Planck’ length ^ Explaining the small magnitude of 'll! is one of the basic 
puzzles of cosmology; here it leads to a condition on physical GFT condensates in order 
to be observationally viable. The resulting picture is very different from the usual one in 
LQG where one thinks of elementary geometric quanta of Planck size, corresponding to the 
lowest non-zero spin j = \l‘l in LQG [46]. In the setting of GFT condensates, (4.8) implies 
that if the average size of the patches is of order k, their relative fluctuations resulting in 
inhomogeneities are at least of order one, so that the universe is very inhomogeneous. 


5 Example: LQC improved dynamics and more 


We have seen in section 3 how global observables, obtained from expectation values of one- 
body operators on the GFT Fock space, can be used to statistically distinguish between the 
homogeneous mode and inhomogeneities. Dynamical equations for such global observables 
can then be interpreted in cosmological terms. Here we give one example of this in a simple 
setting that has been studied before [23-25], but is now reinterpreted. 

The quantum dynamics of a given GFT can be expressed in terms of Schwinger-Dyson 
equations, 




0[ip,ip] 




= 0 


(5.1) 


where 0\(p,ip\ is a functional of the GFT field ip and its complex conjugate ip and S is the 
action for the GFT model. Choosing O = (p{gi) and passing to the canonical operator 


^To make this argument quantitative, one needs more insight on the value of k in the microscopic theory. 
In the standard view of LQG [9-11] (see (3.4) and below), k is given in terms of the low-energy Newton’s 
constant; if k is affected by renormalisation this would lead to a different fundamental ‘Planck’ length. 
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formalism, this becomes 


0H91) 




= 0 


(5.2) 


with normal ordering under which the delta distribution b(plb(p disappears. Integrating 
over we obtain 


K) + 


{dgf(p\gi) 




= 0 


(5.3) 


where K is the quadratic part of the GFT action and V contains all higher order terms. 
The approximation made in [23-25] is now to neglect the second term. This can be an 
exact result for certain states and choices of V [23] , or correspond to a weak-coupling limit 
of the GFT. If we then choose 


^ = f ^\ 9 l) fX] ^{91) 


(5.4) 


which contains a nontrivial propagator as motivated by studies of GFT renormalisation 
[47, 48], comparing with (3.11) and (3.16) we obtain 


'^{ai) - K^g.^{N) = 0. (5.5) 

I 

Focussing on isotropic universes for simplicity, as in section 3, we then identify 



and, assuming all (d/) are equal, the dynamical equation (5.5) reduces to the relation 



between the cosmological variables a, if: and N, where the approximation comes from 
viewing the sum over i as approximating a continuum integral over space. 

Assuming that the universe is exactly homogeneous, (5.7) would imply that 




(5.8) 


which means that for a dynamical state the total physical volume a^Vo is proportional to 
the number of patches N, so that each patch has fixed volume (K/i/2)^/^. 

If (5.8) holds, the ‘average group element’ (II/)/A^ (see (3.15)) for GFT condensates is 
identified with the cosmological observable 


KM 

N 



sin 



(5.9) 


The relation (5.8) reproduces, in GFT condensate cosmology, the improved dynamics pre¬ 
scription of LQG [32] in which holonomy corrections take the form (5.9). In order to obtain 
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(5.9) for GFT condensates, we assumed the form of the kinetic term (5.4), the existence of 
a weak-coupling limit, and consistency of the Poisson brackets at low curvature that led to 
(3.18). Our argument can be contrasted with results in LQC, e.g. [49] where the improved 
dynamics scheme was shown to be the only LQC prescription leading, in the (semiclassical) 
effective dynamics, to a universal ‘quantum gravity scale’ bounding geometric quantities 
that is independent of initial conditions. The improved dynamics scheme has also been 
derived from GFT condensate dynamics in [26] , using a WKB approximation and the cou¬ 
pling to matter. The two results, obtained using different methods, show different ways of 
how LQC holonomy corrections can emerge from the more fundamental GFT framework. 

For the more general case of inhomogeneous universes, (5.7) becomes 

~ (^) ' N (l + ^ J d^x V’(f)^^ (5.10) 

and hence the total volume is decreased when inhomogeneities are present. 

At the end of section 4 we have argued that a nearly homogeneous universe is only 
possible if the length scale associated to the distribution of ‘patches’ is large compared to 
the ‘Planck’ scale. This statement together with (5.10) would imply that we require S> 1 
in order to have, within the approximations we are using here, dynamical condensates that 
are observationally viable, i.e. satisfy 1. This is an example of how input from 

observation can be translated into constraints on a class of GFT models, as was already 
advocated in previous works such as [25]. From the perspective of fundamental quantum 
gravity, (5.4) with fj? ':$> 1 might be viewed as a limiting case in which the Laplace-Beltrami 
operator is a ‘small perturbation’ to the trivial kinetic term on which existing GFT/spin 
foam models related to Plebanski gravity are based, 6 .( 7 . the class of GFT actions in [36]. 

6 Discussion 

The ‘condensate wavefunction’ appearing in the mean-field approximation is not a ‘wave- 
function of the universe’. In deriving quantum cosmology models from GFT condensate 
dynamics, one should use an interpretation that takes the physical interpretation of this 
mean field into account. As we have argued, one way of clarifying this interpretation is 
to rewrite the hydrodynamic approximation of GFT condensates as a classical statistical 
distribution over the space of geometries. This statistical description then already cap¬ 
tures inhomogeneities, and hence this approximate classical picture for GFT condensates 
generally corresponds to an inhomogeneous universe, even though the underlying quantum 
state is homogeneous in a precise sense, being given by a coherent state of many quanta. 

One main future task will be to find explicit solutions to the GFT dynamics, in par¬ 
ticular ones that describe nearly homogeneous universes, in order to be able to extract 
clear predictions from GFT condensate cosmology. It may also be necessary to go beyond 
the mean-field approximation, given that the information about inhomogeneities we can 
extract is rather limited, corresponding to integrals f d^x for different n. 

An important role should be played by solutions that are similar to solitons that appear 
in the physics of Bose-Einstein condensates [31]. The existence of such solutions depends 
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crucially on the nonlinearity of the Gross-Pitaevskii equation. Similarly, the interaction 
term in GFT is the main ingredient for creating Ad spacetime structure out of Sd building 
blocks [10, 12, 13], which also suggests the need for solutions to the full, nonlinear equa¬ 
tions. The analogue of the Gross-Pitaevskii equation for GFT condensates, in the simplest 
approximation where the condensate is defined by (2.6), is just the classical GFT equation 
of motion. Solutions to these classical equations, for GFT models of interest in quantum 
gravity, have been the focus of previous work [50-52] which mainly aimed at exploring the 
dynamics of perturbations around a non-trivial background. One might try to choose one 
of these solutions as a candidate mean field cr{gi) for cosmology. However, the solutions 
used in [50-52] are sharply peaked on ‘flat’ geometries®. This is natural in the context of 
[50, 52] of models for 3d gravity where solutions to the classical theory, defined by a ‘BF 
action’ f tr(E A ^[H]), correspond to flat connections, with conjugate triad 1-form E only 
determined up to transformations E E + d^?? which leave the action invariant [54] . In Ad 
gravity, however, one requires a distribution with finite spread both in terms of curvature 
and metric variables (where the spread in the metric, in order to match observation, must 
be very small), which corresponds to a very different type of classical solution. 

In section 5 we gave one example of how to extract information about the GFT dy¬ 
namics already from the linearised equations. Using the cosmological interpretation of 
expectation values of GFT Fock space operators given in section 3, we have reproduced 
the improved dynamics prescription of LQC holonomy corrections from a consistency rela¬ 
tion between the total volume and the average particle number which means that a® oc N. 
The constant of proportionality depends on the GFT coupling constant fj?] assuming that 
the scale associated to the ‘patches’ in the condensate is ‘large’, aiVo/N)^^^ S> as seems 
required for a nearly homogeneous universe, only models with /r S> 1 would be viable. 

Beyond the simplest condensates we have considered, one might expect ‘large’ patches 
to arise as an effective description at a mesoscopic scale, after coarse graining of fundamen¬ 
tal GFT quanta. The generalised condensates of [29] which involve such a coarse graining 
(see footnote 6) might hence be better suited for describing a realistic cosmology. 

The calculations of section 5 are a special case of the general discussion of [25] where, 
using the variables (3.15), a generalised Friedmann equation was obtained that depended 
on the ‘atomic number’ N which has no classical analogue. Such equations can only be 
interpreted physically when the ‘equation of state’ N = N(a) is known. Here, under the 
assumption of exact homogeneity, we find the relation (5.8) which also depends both on 
the cosmological variable a and the particle number N but simply fixes a in terms of N 
(or conversely), and hence provides the relation N = N(a) but no further information 
about the dynamics which has to come from elsewhere. In a complete derivation of cos¬ 
mological dynamics from GFT, even in the simplest isotropic, homogeneous case, at least 
two independent relations are required to provide these two separate ingredients (similar 
to standard cosmology where the Friedmann equations are supplemented by an equation 
of state). The first step towards this, a main focus of current work [55], is to discard the 

®See also the closely related work in [53] using coherent states, where (again in 3d) imposing the dynamics 
suggests that states should be sharply peaked on flatness, with large spread in the triad. 
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approximation in which interactions are neglected and to instead find (approximate) solu¬ 
tions to the full nonlinear GFT equations. This should also shed light on the dynamics of 
inhomogeneities, which ultimately have to be matched with observational constraints. 
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